一种求解平面接触应力问题的光弹性—数值组合方法

本文提出一种求解平面接触应力最大值的光弹性—数值计算组合方法,仅从光弹性等差线上采集实验数据,然后用牛顿—拉斐逊技术与最小二乘法相结合计算出结果。文中推导了控制方程,阐述了方法,研究了实验数据取值区域、取值个数以及接触长度测试误差等因素的影响。用两类平面接触问题作为应用实例,结果符合工程要求的精度。     

A fourth‐order Hermitian box‐scheme with fast solver for the Poisson problem in a cube

We present a fourth‐order Hermitian box‐scheme (HB‐scheme) for the Poisson problem in a cube. A single‐nonstaggered regular grid is used supporting the discrete unknowns u and . The scheme is fourth‐order accurate for u and in norm. The fast numerical resolution uses a matrix capacitance method, resulting in a computational complexity of . Numerical results are reported on several examples including nonseparable problems. The present scheme is the extension to the three‐dimensional case of the HB‐scheme presented in Abbas and Croisille [J Sci Comp 49 (2011), 239–267]. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 609–629, 2015     

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori–Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.     

二项线性随机效应模型的拟蒙特卡罗估计

讨论了如何运用拟蒙特卡罗方法对二项线性随机效应模型进行参数估计.首先写出观测数据的边缘对数似然函数,然后用拟蒙特卡罗方法将函数中的积分写成求和的形式,接着利用Newton-Raphson算法计算参数的极大似然估计.以一组种子数据为例,说明该方法是简单可行的.     

Un contre-exemple sur les families de courbes gauches

In this paper we give an example of two composants X and X,_o of the Hilbert scheme of space curves (i.e. components of subschemes with constant cohomology) satisfying the condition of semi-continuity (the cohomology of curves in X is less than the cohomology of curves in X,_o) but not the property X ₀?[Xbar]≠Ø.     

变系数单指标模型的B-样条估计

变系数单指标模型结合和单指标和变系数模型的优点,在许多领域中有着重要的应用。本文我们基于B-样条逼近,提出了两种估计方法：第一种方法是利用Newton-Raphson迭代方法同时获得参数和非参数部分的估计;第二个方法是剖面方法获得了相应的估计。当模型中有许多参数时,我们建议使用第二种估计方法,而当参数个数较少时采用第一种方法更方便。两个模拟例子用来验证本文提出的估计方法。     

Improvements of the Maximum Pseudo-Likelihood Estimators in Various Spatial Statistical Models

Abstract

Maximum pseudo-likelihood estimation has hitherto been viewed as a practical but flawed alternative to maximum likelihood estimation, necessary because the maximum likelihood estimator is too hard to compute, but flawed because of its inefficiency when the spatial interactions are strong. We demonstrate that a single Newton-Raphson step starting from the maximum pseudo-likelihood estimator produces an estimator which is close to the maximum likelihood estimator in terms of its actual value, attained likelihood, and efficiency, even in the presence of strong interactions. This hybrid technique greatly increases the practical applicability of pseudo-likelihood-based estimation. Additionally, in the case of the spatial point processes, we propose a proper maximum pseudo-likelihood estimator which is different from the conventional one. The proper maximum pseudo-likelihood estimator clearly shows better performance than the conventional one does when the spatial interactions are strong.     

基于不光滑边界的变系数抛物型方程的高精度紧格式

本文讨论基于不光滑边界的变系数抛物型方程求解的高精度紧格式.首先构造一般变系数抛物型方程的高精度紧格式,并在理论上证明格式具有空间方向四阶精度.然后针对非光滑边界条件,引入局部网格加密技巧在奇异点附近进行不均匀的网格加密.数值实验以期权定价中Black-Scholes偏微分方程的求解为例,验证高精度紧格式用于光滑边界条件的微分方程离散可以达到四阶精度.对于处理非光滑边界条件,网格局部加密技巧能有效的提高数值解精度,使得高精度紧格式用于定价欧式期权可以接近四阶精度.     

牛顿拉夫逊算法在ERT反演问题中的应用研究

本文研究了电阻率反演成像(ERT)中的牛顿拉夫逊基础算法及改进问题.利用最小二乘法和Tikhonov正则化等方法将反演算法予以优化,获得了与实验样本结构吻合的碳纤维复合层的电阻率分布图像,推广了牛顿拉夫逊算法的数理反演模型.